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Abstract. 

Dark matter annihilations in the Galactic halo inject relativistic electrons and positrons 
which in turn generate a synchrotron radiation when interacting with the galactic magnetic 
field. We calculate the synchrotron flux for various dark matter annihilation channels, masses, 
and astrophysical assumptions in the low-frequency range and compare our results with radio 
surveys from 22 MHz to 1420 MHz. We find that current observations are able to constrain 
particle dark matter with "thermal" annihilation cross-sections, i.e. (cry) = 3x 10~ 26 cm 3 s _1 , 
and masses Mom ^ 10 GeV. We discuss the dependence of these bounds on the astrophysical 
assumptions, namely galactic dark matter distribution, cosmic rays propagation parameters, 
and structure of the galactic magnetic field. Prospects for detection in future radio surveys 
are outlined. 
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1 Introduction 

The presence in the Universe of a sizable amount of Dark Matter (DM) has puzzled our 
comprehension about both fundamental physics and the cosmos. One of the most natural 
solutions to this problem is that DM is formed by relic particles in the form of Weakly Inter- 
active Massive Particles (WIMP). WIMPs can in fact easily thermalize in the early Universe 
and non-relativistically decouple from the plasma with the right abundance when they have 
a thermally-averaged annihilation cross-sections of the order of (erf) = 3 x 1CP 26 cm 3 s _1 . 
Dynamically they are then cold enough to allow structure formation as required by the cold 
dark matter paradigm. 

Several ways have been proposed on how to search for WIMP dark matter. Direct 
detection is the most straight technique. Recent results from DAMA [1, 2], CoGeNT [3] 
and CRESST-II [4] all point toward a possible signal induced by a DM with mass in the 10 
GeV range. A different approach is to search for WIMP annihilations or decays that could 
take place in many astrophysical environments like the galactic halo, dwarf galaxies, cluster 
of galaxies, the Sun and the Earth or even during the early dark ages of the Universe. An 
interesting result comes from the FERMI-LAT Collaboration, that recently presented sig- 
nificant bounds on dark matter annihilations into gamma-rays occurring in satellite galaxies 
of the Milky Way [5]. From this analysis, WIMP candidates with "thermal" annihilation 
cross-sections (i.e. (crv) = 3 x 10 -26 cm 3 s -1 ,) and inducing a flux of 7-rays from 7r°-decay 
(i.e., annihilating into quarks or r-leptons) are strongly constrained for masses below few 
tens of GeV [5, 6]. On the other hand, WIMP candidates annihilating into light leptons are 
much less strongly constrained. 
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In this work, we focus our attention on the Galactic synchrotron emission from DM 
annihilations and we concentrate on observations in the frequency range from 22 to 1420 
MHz. We extend previous analyses [7-11] to cover the very low-frequency range, specifically 
we will consider radio emission at 22 and 45 MHz. Such low radio frequencies have not been 
addressed in the past in connection with DM studies, and we will show that they may be 
quite adequate in the search for DM, especially for low-mass WIMPs which are currently 
under deep scrutiny due to the direct detection results quoted above. Indeed, for a magnetic 
field of jiG strength, as typical in the Galaxy, synchrotron emission below GHz-frequency is 
generated by electrons and positrons with energies well below 10 GeV. In particular, we will 
show that WIMP candidates annihilating into light leptons with a "thermal" annihilation 
cross-sections can be strongly constrained for masses Mdm % 10 GeV. 

Radio emission arising from DM annihilation have instead been addressed in Refs. [7— 
11] by using higher frequency surveys and by studying heavier dark matter candidates (more 
suitable to explain the PAMELA positron excess [12]), with a particular focus on the Galactic 
Center by Refs. [13-18], or in order to reproduce the so called "WMAP Haze" in Refs. [19- 
22]. DM searches with radio observations can also be pursued with the cosmological diffuse 
radio background, using both intensity [23, 24] and anisotropies informations [25, 26]. 

The outline of this work is the following. In section 2, we describe general aspects about 
galactic DM and we discuss our approach to the dark matter profiles used in our calculation. 
In section 3, we review electrons and positron propagation in the galactic environment and 
discuss the main sources of uncertainty which will then represent a key element of discussion 
for the predicted radio fluxes. 

In section 4, we discuss the synchrotron emission from a galactic population of electrons 
and describe models for the Galactic Magnetic Field (GMF) distribution. Sections 5 and 6 are 
devoted to the results on the radio emission from DM annihilations and to comparison with 
experimental data. We specifically adopt available temperature skymaps at 5 frequencies: 
22, 45, 408, 820, 1420 MHz [27-31]. From this comparison we then derive bounds on the 
DM particle properties (annihilation cross sections and channels, DM mass). We discuss 
the impact of the uncertainties on the galactic halo modeling, electron propagation, size and 
shape of the GMF. 

In Section 7, we analyze possible traces of WIMP annihilation at 22 and 45 MHz after 
subtraction of a template based on the Haslam map [29] for the astrophysical component. In 
Section 8 we summarize our conclusions. 



2 Dark Matter distribution 

N-body simulations are currently the best tool to study the formation and evolution of 
cosmic structures. Galactic DM halos obtained by these simulations contain a smooth halo 
component together with a population of clumps and other possible unvirialized structures, 
such as streams. 

The radial dark matter distribution of a Milky Way-like DM halo obtained by the Via 
Lactea II (VLII) simulation [32] is well captured by a Navarro-Frenk- White (NFW) density 
profile [33]: 

p NFW (r) = /' x2 (2.1) 

where r is the distance from the center of the DM halo and the values of the parameters 
p s and r s are r s = 21kpc and p s = 0.31 GeV cm~ 3 [34]. The local DM density is then 
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pi = 0.43 GeV cm -3 . The results of the Aquarius simulation are instead better fitted by a 
slightly shallower Einasto profile [35, 36]. In the determination of the NFW parameters we 
closely follow the analysis of Ref. [34], where a detailed study of the dark matter density 
profiles arising from numerical simulations has been performed, and where the relevant profile 
parameters have been determined. With the above values, the halo mass within the virial 
radius is correctly reproduced [34] . 

It is however unclear whether these profiles properly describe the DM distribution in the 
innermost regions of the halo, which are not accessible to numerical simulations due to their 
finite resolutions. Moreover, the feedback of baryons on the DM distribution, not considered 
in the simulations cited above, can be relevant in the inner regions of the Milky-Way. In 
particular, the condensation of gas and stars can produce an adiabatic contraction of the DM 
halo, increasing the central DM density [37, 38]. This effect has been confirmed by recent 
hydrodynamical simulations which include the presence of gas and stars (see Ref. [39] and 
refereces therein). Also the effect of the central super massive black hole on the surrounding 
DM distribution is uncertain since it may either induce an adiabatic contraction of the DM 
halo or destroy the central DM cusp during its spirraling at the center of the halo. The central 
DM cusp can also be erased by hierarchical mergers occurred to the Mily Way halo during its 
recent evolution [40-42] . Current microlesing and dynamical observations are consistent with 
the NFW and Einasto profiles found in numerical simulations and can be used to constraint 
combinations of local DM densities and inner slopes of the DM profiles, ruling out extremely 
cuspy DM distributions [43]. Nevertheless, cored profiles, despite disfavored by simulations, 
are observationally viable and in some cases preferred over cuspy ones, like for observations 
of dwarf spirals and low-surface brightness galaxies [44] . 

In our analysis, in order to bracket the uncertainty arising from the shape of the DM 
profile, we consider the NFW profile discussed above and an isothermal profile, the latter 
used a conservative case of a cored DM distribution: 

P 18 °(r) = po (2.2) 

with a core radius r a = 5 kpc. The normalization po is chosen such that the same dark 
matter density of the NFW profile of Eq. (2.1) occurs at the Sun position, which we take at 
r© = 8 kpc. Therefore po = 1-53 GeV cm~ 3 . 

In addition to the smooth density distribution discussed above, the Milky-Way may 
host a population of DM clumps which would give rise to an additional contribution to the 
electron/positron flux induced by DM annihilations. Unfortunately, the properties of this 
population of subhalos, which are inferred from simulations, are affected by many uncer- 
tainties. Since our study is meant to derive bounds on DM annihilation in the Galaxy, we 
choose not to include substructures in our analysis. In this case, the constraints we derive 
in Sec. 5 are conservative, since the effect of the presence of clumps is to boost up the radio 
fluxes produced by the smooth halo. Nevertheless, we have made an estimate of a potential 
effect of clumps on our results, and we found that their inclusion would not not affect the 
main results of this paper, i.e. the bounds on the DM annihilation cross section. To perform 
this check we have considered a population of clumps modeled following the results of the 
VLII simulation [45-47] and using the prescriptions of Ref. [34]. We found that in a wide 
region (~ 30 degrees) around the galactic center, the average radio flux from the population 
of clumps can be safely neglected as compared to the radio flux produced by the smooth 
halo. This would be even more so for the Aquarius simulation [34] , and occurs because 
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the spatial distribution of the clumps in the numerical simulations occur to be anti-biased: 
the majority of the clumps lie at large radii and their number density is depleted at small 
galactocentric distances. Since the bounds on the dark matter annihilation cross section we 
will derive in Section 6 will mostly arise from the inner parts of the Galaxy, we can conclude 
that for our purposes the average contribution of the subhalos can be safely neglected in the 
present study. The relevant and interesting prospects of detecting individual clumps with 
radio observations is left for a future study. 

3 Cosmic rays electrons 

The propagation of cosmic rays (CR) electrons and positrons is modeled following a semi- 
analytical approach [48]. Cosmic rays are assumed to be confined by the Galactic magnetic 
field inside a propagation zone which is described by a cylinder centered at the galactic 
center. Its radius along the galactic plane R g = 20 kpc approximately matches the radius 
of the galactic disk. The value of the vertical half-thickness (L z ), which is instead more 
uncertain, is estimated from nuclear cosmic rays observations and takes values from 1 to 
15 kpc [48, 49]. Outside the propagation zone, cosmic rays are no longer confined by the 
magnetic field so their density is expected to rapidly drop. The CR propagation is described 
by the transport equation for the number density of electrons per unit of energy (ip): 

- K E s V 2 iP - A = g ( X) E ) , (3.1) 

where Kq and 5 are parameters related to spatial diffusion, b(E) is the energy loss term, and 
q(x, E) is the source term, that in our case will correspond to DM annihilations. 

The energy loss term for electrons takes into account the synchrotron and inverse Comp- 
ton (IC) scattering losses from interactions with magnetic and interstellar radiation fields 
(ISRF): 

b(E) = b 1SRF (E) + b synch (E, B) , (3.2) 

where the term related to the ISRF, &isrf> is modeled as in Ref. [50]. We take under consid- 
eration a specific modeling of the various ISRFs and we explicitly include the Klein-Nishina 
corrections to Compton scattering [50]. Nevertheless, these corrections do not produce such 
important changes in the final result specially for energies lower than 10 GeV, where Comp- 
ton scattering is in the Thomson regime. As a reference value, for electron energies below a 
few GeV (which are relevant for light dark matter), &isrf — 0.8 x 10 -15 (i£/GeV) 2 GeV/s. 
The energy losses from synchrotron emission depend on the galactic magnetic field intensity 
B as [51]: 

4e 4 -B 2 

b sy nch(E,B) = gm 4 c7 E , (3.3) 

where c is the speed of light, e and m denote the electron's charge and mass, respectively. 
For a magnetic field of 0(6fiG), the ratio of the synchrotron and IC losses is: 

frsynch/klSRF ^ 1/2.5 , 

which indicates that 6 S ynch does not dominate the full energy loss term. For this reasons, for 
conservative values of B(< 20 fj,G) the magnetic field can be safety be assumed to be constant 
inside the propagation zone [50, 52]. This assumption is well justified for what concern the 
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Model 


propaga 

L z [kpc] 


:ion param 

^ r k P c2 

K Myr 


Leters 

5 


MIN 


1 


0.0016 


0.85 


MED 


4 


0.0112 


0.70 


MAX 


15 


0.0765 


0.46 



Table 1. Benchmark propagation models compatible with B/C and p/p observations [48, 49, 61]. 

CR propagation while the morphology of the associated synchrotron emission depends on 
the details of the galactic spatial structure. 

In the galactic disk (~ 200 pc thick) additional energy losses, like bremsstrahlung and 
gas ionization, might become relevant, in particular at low electrons energies [53]. While 
these effects might be important for CR sources localized on the disk, they are expected to 
play a minor role in the case of more extended sources, like a DM halo, significantly reducing 
the synchrotron flux only at latitudes \b\ < 1°. 

The parameters of the propagation model are constrained by CR data, notably the 
Boron over Carbon ratio (B/C) and radioactive isotopes [48, 49, 54]. Let us mention that 
the resulting viable parameter space is also consistent with the secondary positrons data [53] . 
In this study, we use the following three propagation benchmarks: MIN, MED, and MAX (Ta- 
ble 1). These models are fully compatible with B/C and antiproton-proton ratio (p/p) ob- 
servations, describing the extremes behaviors (MIN and MAX) and best fit (MED) of these 
observables [55, 56]. 

Recent efforts to further constraint the space of the propagation parameters have focused 
to "alternative" observables, like synchrotron total intensity [57, 58] and anisotropics [59], 
and gamma rays [60]. 

3.1 Galactic electron population from dark matter 

In addition to standard CR astrophysical sources, dark matter annihilations in our Galaxy 
could provide an exotic source of CR electrons and positrons. The resulting electrons number 
density is computed from the transport equation Eq.(3.1) with the following source term: 



g(x,£7) = 




where (crv) is the annihilation cross section, p(x) is the dark matter density and Mdm is the 
dark matter mass. The function dn/dE is the energy spectrum of electrons/positrons per 
single annihilation and depends on the particular DM annihilation channel [62, 63]. 

The transport equation can be solved using the formalism of Refs. [62, 63], which is 
based on a decomposition of the source term in Fourier and Bessel modes, exploiting the axial 
symmetry of the propagation zone (for more details see Ref. [62]). This method, initially 
proposed to calculate the electrons/positrons fluxes at the Earth position, can be extended to 
compute the electron density anywhere inside the propagation zone. The spatial information 
of the source term is encoded in the so-called halo function, / x , which for a particular position 
x depends only on one variable Ad: 

L(E,E S ) = L(X D (E,E S )), (3.5) 
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Injection energy [GeV] 



Figure 1. Diffusion iength Ad for propagated electron of 0.3 GeV versus injection energy for 
MIN, MED, MAX propagation models (Table 1) and magnetic field intensities of 6 /j,G and 10 /xG. 
Depending of the propagation model, electrons can easily cover regions larger than kpc 3 . 



where Ad is the diffusion length 



\UE,E S ) = 4 de-^-, (3.6) 

The diffusion length depends on the propagation model and energy losses and it can be under- 
stood as the scale length that can be reached by electrons with energy E that were produced 
with energy E s . For propagation models compatible with B/C observations, electrons can 
cover distances larger than few kpc's becoming mildly insensitive to small-scale variation of 
energy losses (Fig.l) which can be therefore approximated to be spatially constant. 
The source term in Eq.(3.4) can be factorized as: 

q(x,E) = f( X )xg(E). (3.7) 

The electron number density is obtained as: 

M*> E) = — dE s I X (E, E s ) g(E s ) , (3.8) 

where the halo function at position x(= (r, z)) is give by: 



R, 



Jx(A D ) = ^ ^ C in Jo (-^r ) <j> n {z) exp ( ™ ) ' ( 3 - 9 ) 

i=l n=0 ^ 9 



The coefficients Ci n are obtained through the decomposition of /(x) in Fourier and Bessel 
modes. Jo(x) is the first kind bessel function, ao,i is the i-th zero of Jo and (fi n are the 
Fourier eigenfunctions: 

a ( z \ = { sin ( fcnZ ) n even f3 ioj 

n | cos (k n z) n odd 
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GMF Model 
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L m [kpc] 


leters 

Rm [kpc] 


I 


5L Z 


6R g 


II 


L z 


Rg 


III 


1 


Rg 


IV 


constant 



Table 2. Parameters of the galactic magnetic field models considered in this study. We assume that 
the magnetic field at Sun position (r Q = 8kpc) has a value of 6 /J.G. 



with k n = mr/2L z . Finally, Ki n is given by: 

i^ = (|£) 2 +e (3.H) 

For well behaved functions /(x), the halo function is well defined for = 0, where it 
converges to the value: 

4(A D -> 0) -> /(x) (3.12) 

To ensure the numerical convergence everywhere in the propagation zone and for diffusion 
lengths larger than 10 pc, we needed to calculate up to 10 7 coefficients Cj n . 



4 Synchrotron radiation 

Electrons propagating inside the Galaxy interact with the galactic magnetic field (GMF) 
producing synchrotron radiation [51]. For magnetic field intensity of 0{^G) [64], like in the 
case of our galaxy and for electrons/positrons of GeV-TeV energies, the synchrotron emission 
falls in the MHZ-GHz range [65], i.e. in the radio band. 

The synchrotron flux observed at the Earth results from the integration along the line 
of sight of the synchrotron emissivity j u : 

j„(x) = J dEn e (E,x)^(x) (4.1) 

where n e (E,x) is the energy spectrum of the electronic number density and dw/dv denotes 
the emission power, which depends on the magnetic field strength -B(x) (see Appendix B for 
details). In our case, the electron spectrum is obtained by solving the diffusion equation for 
the population of electrons and positrons produced by DM annihilations, as described in the 
previous section. 

The structure of the GMF is still not well understood. Mainly, it is composed by a 
regular and a turbulent part, the last one responsible for the diffusive behaviour of the cosmic 
rays. The main methods used to estimate the GMF are Faraday rotation measurements of 
pulsars (more sensitive to the regular part) and diffuse synchrotron emission from the non- 
thermal electron population (sensitive to both components). 

In the literature, various parametrization of the GMF distributions have been pro- 
posed [66-70]. Here, we consider a magnetic field with a cylindrical symmetry and an expo- 
nential dependence on the distance on the galactic center r: 
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B(r,z) = £ exp (- L ^) /(*) ■ (4-2) 

The regular component of the GMF typically has -Bo ~ 2/iG and a scale radius R m ~ 8.5 kpc. 
Its vertical dependence, encoded in f(z), is more uncertain, due to the low number of pulsars 
at high latitudes. Uncertainties on the turbulent component of the GMF arise from the 
fact that the synchrotron emission, used to estimate it, depends both on the magnetic field 
intensity and the electrons density. This degeneracy can be partially broken using local 
cosmic rays observations. 

An interesting link between the propagation model and the magnetic field is given by the 
expression of the diffusion term of the transport equation Eq. (3.1) for isotropic turbulences 
in the quasi-linear approximations [71, 72]: 

D(E)=K E s otr* I §P±, (4.3) 

"turb 

where r g = TZ/B meSin is the Larmor radius, which depends on the rigidity (TZ), B mea _ n and 
-E>turb are the mean and the turbulent components of the GMF. This allows to connect the 
spatial extension of the GMF with the diffusion coefficient. 

In the following we model both the regular and turbulent components of the magnetic 
field with a double exponential: 

B(r,z) = B exp {--^ - ^] • (4.4) 

If we assume for the scale parameters R m and L m the following relations: 

L m = 5L Z , R m = 5R g (4.5) 

where 5 is one of the diffusion parameters (Eq. 3.1), then the resulting diffusion term (Eq. 4.3) 
will grow exponentially with scales R g and L z .This exponential dependence of the diffusion 
term on the cylindric coordinates is commonly adopted in numerical codes like Galprop [73- 
75] and Dragon [76, 77]. Inside the propagation zone, this parametrization is almost equiv- 
alent to a constant diffusion term. In addition to this parametrization we have considered a 
double exponential model with the same scales as those of the propagation zone and the sim- 
plistic case of a constant magnetic field. As a final additional choice, we have fixed the vertical 
scale of the magnetic field L m to 1 kpc, which is much lower than the height of the diffusion 
zone in the case of the MED and MAX propagation models. These four parameterizations are 
summarized in Table 2. The strength of the total magnetic field (turbulent plus regular) at 
the Sun position is set to B = 6 /j,G, which is well consistent with the measurements [66-70]. 



5 Comparison with observations 

Using the formalism explained in the previous sections we have generated skymaps of the 
synchrotron radiation produced by DM annihilations at different frequencies. The three 
examples shown in Fig. 2 are for a 10 GeV DM particle annihilating into muon pairs and 
for the three propagation models of Table 1. Throughout this Section, for definiteness the 
DM annihilation cross section is set to the thermal value: (o~v) = 3 x 10 -26 cm 3 s . In 
Fig. 2, the NFW density profile and the GMF model I are used. The morphology of the 
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v [MHz] 


Survey 


rms noise [K] 


22 


DRAO [27] 


5000 


45 


Guzman et al.[28] 


3500 


408 


Haslam et al. [29] 


0.8 


820 


Dwingeloo [30] 


1.4 


1420 


Stockert [31] 


0.02 



Table 3. Surveys considered in the analysis. The last column is the rms temperature noise. 

emission strongly depends on the propagation model, DM spatial profile, and the structure 
of the GMF. In particular, the propagation zone height L z has a large impact, with larger 
values of L z leading to more extended emissions. 

Recently, radio observations have been used in combination with cosmic rays measure- 
ments at Earth to constraints the cosmic rays propagation model [57, 58]. Here, we instead 
study the bounds on the dark matter annihilations cross section which can be inferred from 
current data and we investigate the possibility to detect dark matter annihilations with radio 
observations. 

From the different radio surveys available, we select those in Table 3 and plotted in 
Fig. 3. They all refer to observations at < GHz frequencies and have a large sky coverage. In 
Fig. 4 we show the average temperature observed at the galactic poles at different frequencies. 
The data have been averaged over a 10° circle. For illustration, we also plot the temperature 
as a function of the frequency obtained using the software presented in Ref. [78], which 
produces skymaps of the radio emission in a large range of frequencies interpolating between 
available data for the frequencies and the regions of the sky not covered by observations. 

The vertical bands illustrate the average frequency of the synchrotron radiation pro- 
duced by electrons of a given energy and suggest that the the synchrotron emission produced 
by 100 GeV) DM candidates peaks at sub-GHz frequencies. These frequencies are there- 
fore particularly suitable to search for light/intermediate DM masses. Moreover, both WIMP 
models with dominant hadronic/bosonic annihilation final state and leptophilic light WIMP 
models induce a synchrotron spectrum which is softer than the observed galactic one, as 
shown in Fig. 4. One of the main motivations of this paper resides therefore in the fact that 
low radio frequencies have the largest constraining power for such WIMP models. 

The plots in Fig. 5 compare observations and DM emissions at 45 MHz for different DM 
models and astrophysical setups. The observational data correspond to a thin strip crossing 
the galactic center and perpendicular to the galactic plane (|/| < 3°). We notice that the 
emission is more extended for propagation models with larger L z . This means that for the 
MED and MAX models it is possible to search for the DM signal outside the inner galactic 
center region (\b\ > 10°). Inside this region, depicted by the blue band, the astrophysical 
uncertainties on the propagation model, the DM distribution, and the contamination from 
the astrophysical background become more important. As expected, for the cored isothermal 
profile, the emission is strongly reduced at small galactic latitudes and longitudes. The 
scaling of the signal with respect to the magnetic field strength is oc B~ 2 . In general, the 
emission does not dramatically change for different GMF spatial profiles, unless one considers 
a large mismatch between the size of the diffusion box and the scale at which the magnetic 
field becomes suppressed, as it happens for the GMF IV together with the MED and MAX 
propagation models. 

While the observed radio emission is dominated by the astrophysical background, Fig. 5 
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Figure 2. Temperature skymap of the galactic radio emission observed at 45 MHz (upper left) and 
synchrotron skymaps at the same frequency from DM annihilations for the three propagation models 
in table 1: MIN (upper right), MED (lower left), and MAX (lower right). The DM synchrotron 
emission is computed for a 10 GeV DM particle annihilating into muon pairs with an annihilation 
cross section (av) = 3 x 10~ 26 cm 3 s _1 . We adopt the GMF model I and the NFW profile discussed 

in section 2. Color scale corresponds to log 10 ^T/K (j//MHz) 2 ' 5 j and it is common for all maps. 

suggests that DM could substantially contribute to the radio flux, especially close to the 
galactic center region. Certainly a possible DM detection is challenged by the large uncer- 
tainties which affect the determination of the background. In the next section, in order to 
be conservative we just use the present observational data to constrain DM models, without 
attempting any background subtraction, and we comment on the possibility to detect DM 
signals with radio observations in Section 7. 

6 Constraints 

The bound on the dark matter parameter space are conservatively inferred requiring that 
dark matter synchrotron emission does not overshoot the data, without taking into account 
the contribution of the astrophysical background. Our procedure is the following: the sky is 
divided in 232 patches obtained as explained in Appendix A. Then we compute the average 
DM radio emission (Tdm) and average observed radio flux (T b s ) in all the patches which are 
covered by the data (note that most of the surveys do not cover all the sky, see Fig. 3). The 
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Figure 3. Temperature skymaps (from top to bottom and left to right) at 22, 45, 408, 820, and 1420 
MHz. See table 3. Color scale is as in figure 2. 



constraint inferred from each patch is obtained requiring that: 

Tdm < Jobs + 3cr 

where a is the rms temperature noise that we take as in Table 3. We repeat the calculation 
for all the frequencies and then we set a bound on the DM annihilation cross section (cry) by 
taking the most constraining patch. The results are summarized in Figs. 6 and 7. 

The bounds do not dramatically vary for different annihilation channels. We find that 
models with DM masses Mdm ^ 10 GeV and thermal value of the annihilation cross section 
(cry) = 3 x 10~ 26 cm 3 s^ 1 are strongly constrained in the case of the NFW profile. We also 
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log 10 (v/MHz) 

Figure 4. Temperature versus frequency calculated for the galactic poles (6 = ±90°) for and 
bb annihilation channels, Mdm = 10, 100 GeV. The DM profile and the galactic magnetic field are as 
in fig. 2 and the propagation model is the MED one, table 1. The data points are the temperature at 
north and south galactic poles averaged in a 10° circle. Green dashed lines are linked to observations 
and have been obtained with the software developed in [78] (see text for more details). 

notice that the constraints drastically weaken for the isothermal profile, which presents a 
much lower DM density in a large region around the Galactic center. This effect has also 
been shown in Fig. 5. Still, for this profile and for the [i + [i~ and e + e~ annihilation channels, 
thermal values of (av) are excluded for Mdm ^5 4 — 6 GeV. 

Despite the morphology of the emission is quite different for the three propagation 
models, the derived bounds are instead similar (except when cutting away \b\ < 15°), as it 
is shown in Fig. 6, where the case of DM annihilations into muons is considered. This is 
because in the most constraining patches (i.e., low latitudes) the average DM emission is 
similar for the three cases. 

Similar conclusions can be drawn about the GMF model. Indeed, changing the magnetic 
field profiles does not dramatically alter the constraints, since at low latitudes the different 
models considered for the GMF are similar. The example shown in Fig. 6 is for the MED 
model. Similar results are obtained for the MAX setup while in the MIN case bounds are 
weaker, but within a factor of less than 1.5. 

We also study the constraints inferred individually from each survey (upper-right panel 
of Fig. 6 ). As discussed in Sec. 5, the lowest frequencies are more constraining for low 
DM masses (Mdm ^ 10 GeV) while O(GHz) frequencies becomes relevant for heavier DM 
candidates. Let us remark that the constraining power of a single survey also depends on 
the fraction of sky covered and on the sensitivity of the map, as it is shown in Fig. 6 where 
the 820 MHz survey provides worse constraints than the 1420 MHz one. 

As commented in Sec. 3, additional energy losses than those we have considered here 
might become relevant in the galactic plane. We estimate that this effects might reduce our 
predictions on the synchrotron fluxes in the region \b\ < 1°, so their impact on the bounds 
would be rather small, since they are derived considering patches of the sky with significantly 
larger angular sizes. Still, we decide to compute the bounds by cutting a large region around 
the galactic plane by imposing \b\ < 15° (lower-left panel in Fig. 6), which in our case is 
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Figure 5. Temperature versus galactic latitude at 45 MHz. The data corresponds to a thin 
strip |/| < 3°, with / the galactic longitude. Lines are predictions for DM models for / = 0°. Blue 
band corresponds to |6| < 10° and it denotes the directions towards to the Galactic center region. 
Upper-panel shows prediction for GMF models described in table 2, for NFW profile, MED model, 
and channel. Middle-panel shows predictions for different annihilation channels, NFW and 

Isothermal profiles with GMF set to model I and propagation to the MED model. Lower-panel shows 
predictions for the channel with the all three propagation models: MIN, MED, and MAX, and 

with the profiles: NFW and Isothermal. 
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10° 10 1 10 2 10 3 10° 10 1 10 2 10 3 



M DM [GeV] M DM [GeV] 

Figure 6. Upper bounds on the annihilation cross section (av) as a function of the DM mass for the 
muonic channel. Different panels are for different astrophysical setups. Upper-left panel shows the 
constraints inferred with all survey using the whole sky. The constraints for NFW profile are more 
stringent than the Isothermal case. Upper-right panel shows the constraints inferred individually 
from each of the surveys considered in the analysis. In the lower-left panel, the bounds are obtained 
excluding the strip along the galactic plane \b\ < 15°. This means that in the grid of appendix A 
we have excluded three lines of patches around the galactic plane. Constraints inferred for different 
GMF models are presented in the lower-right panel. 

equivalent to remove 3 strips of patches. In general, with this conservative choice, we reduce 
the uncertainties on the propagation related to possible unaccounted astrophysical effects 
occurring on the galactic disk. Moreover, the fact that the radio astrophysical background 
and its uncertainties are maximal in the galactic plane, complicates the searches of DM signals 
in this region. Therefore, by avoiding the region \b\ < 15°, we study the reach of current 
radio surveys in the DM parameter space focusing on "cleaner" targets of observations. 
Interestingly, we find that the bounds are almost unchanged for the MAX model and only 
slightly affected for the MED one. For the MIN model, they are more altered since the 
emission is more concentrated in the Galactic center region. These results are consistent 
with what discussed in the Sec. 5. 

7 Prospects for detection 

So far our goal has been to derive robust and conservative bounds on WIMP parameter space. 
We now analyze the possibility that the data already contain a significant DM contribution, 
and discuss a method to single it out. 

As for the constraints, our focus is on WIMP models inducing a synchrotron emission 
with a spectrum softer than the astrophysical one (see Fig. 4) with the latter mainly coming 
from cosmic-rays and showing a spectrum roughly proportional to v~ 2 ' 5 . As it is clear from 
Fig. 5, a possible DM excess would be maximal in the central region of the Galaxy (within 20- 
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M DM [GeV] M DM [GeV] 

Figure 7. Upper bounds on the annihilation cross section (av) as a function of the DM mass 
for different annihilations channels: e + e~ (Upper-left), t + t~ (Upper-right), bb (Lower-left), and 
W + W~ (Lower-right). In most of the cases, the thermal cross section value, orange solid line, is 
reached at Mdm ~ 10 GeV. 



30 degrees) , and should appear as a spherical feature related to the (approximately) spherical 
shape of the DM profile. This is opposite to what is expected for the CR emission, which 
typically shows a "disky" shape, following from the confinement of CR sources to the stellar 
disc. Therefore, we perform a search where the signal to background ratio is expected to be 
larger, i.e. in maps at low frequencies and in the inner part of the Milky-Way. We choose 
the 45 MHz map [28] where the central region is better sampled, but analogous analysis can 
be done with the 22 MHz map [27] . 

We consider the Haslam et al. map [29] at 408 MHz as a template for galactic syn- 
chrotron emission. This is commonly done also in CMB studies since the Haslam map is the 
radio full-sky map with the best angular resolution and sensitivity, at a frequency where the 
emission is thought to be dominated by synchrotron radiation. We estimate the emission at 
45 MHz in each angular pixel as Tf st = T 408 ,i ■ (45MHz/408MHz) a + T°, where i is the pixel 
index, and the 408 MHz map is smoothed down at the angular resolution of the 45 MHz 
map [28]. a and T° are derived from the best-fit of the observed temperature T obs in the 
map of Ref. [28], by minimizing the x~ s 1 uave function \ 2 = Si(^T st — T° hs ) 2 /a 2 , where a 
is the noise-level reported in Table 3, and in the sum we include all available pixels in the 
map except for the disc (|6| < 5°). We include T° in the fit to account for a possibly different 
spectral index in the extragalactic emission, or for possible experimental issues associated to 
the absolute normalization of the flux. However, it is nearly irrelevant since its best-fit value 
comes out to be at the level of noise. We found a = —2.56, confirming our expectation, and 
suggesting that this search technique can be indeed useful for DM candidates inducing a syn- 
chrotron spectrum softer than ~ v~ 2 ^ '. The residuals (T? st — ^"B^rpobB are shown in Fig. 8 
for the central 30° x 30° box, and in the full-sky. Note that this technique is similar to the one 
employed for the "WMAP haze" [19, 20, 22], namely a possible foreground excess found in 



- 15 - 



the inner Galaxy at microwave frequencies and which can be interpreted in terms of WIMP 
annihilations/decays. Contrary to the possibility investigated in this Section, the WMAP 
haze involves an excess which is harder than the rest of the galactic synchrotron emission. 
Therefore, in that case, the focus is on high frequencies and different DM candidates. 

No evidence for a bright and spherical spot is found in Fig. 8, with the residuals fluctu- 
ating between ±30% in the central region. Therefore we found no hints for a DM signal in 
current data. Similar conclusion can be drawn considering the 22 MHz map [27]. However, 
as clear from Fig. 8, our estimates of the temperature are not matching well with measured 
ones, with large residual fluctuations which are left after subtraction all over the sky. From 
a physical point of view, even assuming that the emission is completely due to synchrotron 
radiation, one can expect a spatial variation in the spectral index (which is assumed to be 
constant in our template technique), associated to variations of galactic magnetic properties 
or energy losses across the Galaxy (which in turn affect cosmic-ray spectrum and emission). 
This should be correlated to the galactic disc, while fluctuations in Fig. 8 are not. Therefore 
they should have a different origin. Most of the brightest spots in the full-sky image can 
be ascribed to known emissions (LoopIII, Cignus region, Cen A, Virgo, North Polar Spur, 
and Gum Nebula to mention the most relevant) other from the diffuse galactic synchrotron 
radiation. However, there is still some large residual left (see, e.g., bottom-right corner in 
the central region map). 

On the other hand, both the derived spectral index and residuals are in fair agreement 
with results found in Ref. [78]. Adopting a principal component analysis, Ref. [78] estimated 
that the first principal component at low radio frequency can be in first approximation 
identified with a synchrotron component with a spectral index a ~ —2.5 and accounting for 
~ 70% of the total emission, while the remaining ~ 30% arises from other components. 

The large residuals can be therefore due either to other physical contributions (i.e., 
free-free emission, spinning dust or thermal dust), although they are expected to be largely 
subdominant, or to possible spurious effects due to poor angular resolution and sensitivity 
(e.g., one can see, even "by eye" from Fig. 3, that the emission in the map at 45 MHz 
is systematically more diffuse than in the 408 MHz map). For these reasons, our test is 
somewhat inconclusive, simply telling us that there is some extra-components which account 
for ~ 20 — 30% of the total emission, but with properties different from the ones expected 
for DM. 

To fully perform this test in the direction of DM searches, one would need more refined 
maps. Future observations with the LOFAR telescope [79] have intriguingly larger possibili- 
ties of detection. Indeed, LOFAR can go down to 15 MHz (thus enhancing the possible signal 
to background ratio for DM candidates discussed here) with unprecedented sensitivity, and 
sample the central region of the Galaxy [80] with a field of view of 20-50 degree and angular 
resolution of 1-30 arcmin, depending on the configuration. This will allow to significantly 
improve the template fitting subtraction and can lead to significantly lower residuals, which 
might help to single out a DM contribution at few-percent level. 

8 Conclusions 

We have argued that radio observations in the low-frequencies regime are particularly suitable 
to search for DM annihilations in the galactic halo. To show that, we have compared available 
radio surveys with the synchrotron emission produced by DM annihilations for different DM 
setups and astrophysical assumptions. We have considered DM masses in the range 1 GeV 
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Residuals 45 MHZ - full Sky Residuals 45 MHz - inner Galaxy 




Figure 8. Residuals R. t = (T° st - T° hs )/T° hs at 45 MHz for the full-sky (left) and in the inner 
30° x 30° of the Galaxy (right). Haslam ct al. map [29] at 408 MHz is used as a template to subtract 
galactic synchrotron emission. Sec text for details. 

- 1 TeV, various annihilation channels which imply different spectral features and, in order 
to bracket uncertainties in the spatial profile, two DM density distributions: the NFW and 
the isothermal profiles. Uncertainties related to the propagation of cosmic rays are taken 
into account considering three benchmarks propagation models, the so-called MIN, MED 
and MAX models. The morphology of the corresponding synchrotron emissions are quite 
different, mostly because these models have rather different sizes of the diffusion box. Indeed, 
larger values of the half-height of the propagation zone L z produce more extended emissions. 
We have then explored different shapes of the GMF and determined their impact on the 
morphology of the synchrotron emission. 

Comparing the DM synchrotron fluxes with low-frequencies radio survey we have com- 
puted upper bounds on the DM annihilation cross section (<rv). Interestingly, we found that 
for a NFW profile a "thermal" cross section is strongly constrained for AfrjM ^$ 10 GeV. The 
constraints we have obtained are similar or stronger than those obtained with other indirect 
DM searches, like gamma-ray observations of dwarf-galaxies [5], galactic center region [81], 
galaxy-clusters [82] and isotropic diffuse background [83, 84]. The bounds mildly depend on 
the propagation model and GMF profile, while they are significantly relaxed for shallow DM 
distributions, like for the isothermal profile. 

In the last part of the paper, we have tried to single out the presence of a possible DM 
component in the present data. We considered the 408 MHz map as a spatial template to 
estimate the galactic synchrotron emission in the low frequency maps. This method allows 
to search for the presence of a DM-induced component softer than the astrophysical galactic 
synchrotron emission, and the DM signal would appear as an excess in the central region with 
an approximately spherical shape. We have concluded that present data does not support 
any evidence for the presence of this additional synchrotron component. However, we expect 
that future radio surveys, in particular with the LOFAR telescope [79, 80], will improve 
coverage, angular resolution, and sensitivity in low frequency radio maps, which may allow 
to disentangle a faint DM contribution. 
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A The grid 

The divide the sky with a grid in spherical angles 9 and 4> which depend on the following 4 
parameters: 



• # ff sets angular scale of poles. 

• rig sets number of #-belts (without counting polar caps). 

• n^ k sets number of ^-patches at the galactic plane (9 = n/2). 

• n^ 1 sets number of ^-patches for the belts closest to the poles. 

In this work, we choose o g = 7.5°, rig = 18, n^ k = 22, and n^ 1 = 6. This specific 
configuration produces a grid where the galactic plane (9 = vr/2) is covered with one belt in 
9. Moreover, it produces two columns centered on <f> = and (j) = it covering the galactic 
center and the galactic anticenter (figure 9). 

The first step is to generate the #-belts. Each one is bounded by 9i and 0i+\, defined 

by: 

for i = 

9 oS + (tt - 29 oS ) ^\ for 1 < i < n e (A.l) 
n e — 1 

7r for i = rig + 1 
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There are in total ng — 1 belts, that cover the region 8 G [0 o s, vr — o g\, and 2 patches covering 
the poles. Each #-belt is divided as well into j' max -+- 1 equally distributed patches in eft: 



^off 



for j = 



6 off + (2vr - 2</> off ) J aY 1 - for 1 < j < f 



(A.2) 



JJLU 



l 



where <^> g = ir/jf 1 ^ ensures a column of 0-patches centered on cj) = 0. 
These parameters also control the value of j'P 1 ^: 

<+(< k -<)i^ I forl<z<^ 



»lM<-<)£-|for J < i < n e 

2 1 



(A.3) 



which linearly interpoles (using division of integer numbers) between n^ 1 and n^ k . 

Then, each skypatch is defined by 4 points in the 9—cj) plane. We follow the convection 

= ($min> $max) x (^min; ^max) j (A-4) 



where k is just the patch index, k is assigned by counting the number of patches going sys- 
tematically from the North to the South pole by covering first each <j!>-patch in every #-belts. 



The Northen polar cap is: 

= (0, 6 oS ) x (0, 2vr) , 
and the patches in the North hemisphere are: 

for 1 < % < n e /2 and < j < jf 3 * - 1. 

The patches in the South hemisphere are 

&k = (0i-i,0i) x ( ( /) ij -,0j iJ - +1 ) 
for n /2 + 2<i< n e /2 and < j < jf 13 * - 1. 

The last patch corresponds to the Southern polar cap: 

Cl k = (vr - 6» ofr ,vr) x (0, 2vr) . 

B Synchrotron emission power 



(A.5) 
(A.6) 



(A.7) 



The synchrotron emission power at frecuency v produced by a relativistic electron with energy 
E [51] is: 

/.■>■ -i ,-, 

v 



dw V3e 3 B ± 

— (v,B ± ) = 2 — F 

dv m e c z \ v r 



(B.l) 
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where 

3eB ± E 2 



47rm^c 5 



/>oo 

and F(x) = x d( K 5/3 (() , (B.2) 

J X 



with -^5/3 is a modified Bessel function. The emission depends only the perperdicular com- 
ponent of the magnetic field (B± = B sin8) with respect to the electron's momentum (9). 

In the case of cosmic rays, the dependence on 6 will be softened due to the diffuse 
(isotropic) propagation of cosmic rays. In this case, we define the function: 

7T 

G( x ) = - [ 2 desm6F(^—) (B.3) 



7T In \Sm 



which is the angular average of the F function. Moreover, the G function can be well fitted 
by: 



G{x) = ax d exp ( - , (B.4) 

where the numerical values of parameters a, b, c, d are: 

a = 1.60883, b = 1.95886, c = 1.13147, d = 0.33839 . (B.5) 
Using this function, the emission power for an isotropic distribution of electrons is: 

dw V3e 3 B 



v 

dv m e c 2 V v, 



, (■{ .. ) • (B-6) 



where v c = 3eB E 2 / (Anm^c 5 ) and B is the total magnetic field intensity. 
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